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1 Introduction. 



In the past few years huge progresses have been made in the study of of nuclei 
at zero temperature in an ab-initial fashion, using both phenomenological effec- 
tive nucleon-nucleon interactions or effective interaction derived from more fun- 
damental bare nucleon-nucleon interactions (see for example refs. [1],[2] and 
references in there for recent reviews). Unfortunately for nuclei at finite tempera- 
ture no such advances have been reported so far. Most of the studies reported so 
far in the literature have been made for simplified phenomenological interactions. 
In the Shell Model Monte Carlo approach usually interactions with good sign 
statistics are used (see for example ref. [3], [4] and references in there). Or with 
approximate theoretical schemes. For example the static path approximation, al- 
though formulated for general interactions has been applied only for quadrupole- 
quadrupole or pairing+quadrupole phenomenological interactions (ref. [5] but not 
to the more complicated realistic effective interactions such as FPD6 or GXPFl A 
or similar, or to effective interaction derived from the more fundamental bare 
nucleon-nucleon interactions. This is not very surprising since even at the level of 
the static path approximation the inclusion of all multipole parts of the interaction 
leads to statistical density operators that are not hermitian and therefore the eval- 
uation of the partition function can be rather difficult, if possible. Surprisingly, 
even the Hartree-Fock-Bogoliubov with the restoration of all quantum numbers 
(cf. ref. [6] for a recent reformulation) has never been applied at finite tempera- 
ture with effective interactions (whether realistic, that is fitted to the experimental 
data, or derived from nucleon-nucleon interactions with renormalization meth- 
ods). The purpose of this paper is to analyze some of the difficulties associated 
with the thermal HFB approximation (THFB) and to propose a simple method to 
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work directly in a canonical ensemble (that is all traces are taken in the subspace 
of the full Hilbert space that has a specified total number of particles). The use of 
quantum-number projected statistical density operators is very relevant since they 
smooth out sharp phase transitions (ref.[7]). This paper is organized as follows. 
In section 2, we briefly review the projected HFB approximation at finite tem- 
perature. In section 3 we summarize the main properties of the statistical density 
operators, discuss the methods needed to evaluate traces in a canonical ensemble 
as well as gradients (with respect to the variational parameters) and expectation 
values. In section 4 we present some concluding remarks. 



2 Recap of the projected THFB approximation. 

There are two methods to obtain the THFB with restoration of exact quantum 
numbers. One method consists in the minimization of the free-energy defined as 

F^<H> -TS (1) 

where H is the many-body Hamiltonian, and 

<H>= Tr{wH) (2) 

with w being the statistical operator 

w^Pe-^^op/Z (3) 

T is the temperature T — 1//3, Hq is a trial quadratic Hamiltonian in the creation 
and annihilation operators that breaks all symmetries, including conservation of 
particle number, P is the projector that restores good quantum numbers, and 

Z = TriPe-^"") (4) 
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is the projected partition function for this trial Hamiltonian. S is the entropy 

defined as 

S — —Tr{wln{w)) (5) 

The trace operation is taken over the full Hilbert space. The minimization is car- 
ried out on the trial Hamiltonian ^o- The difficulty with this approach is the pres- 
ence of the projector P in the logarithm. It can be remedied by replacing, as an ap- 
proximation (cf. ref.[6]), the entropy with /3Tr{ti!Ho)-\-\n Z. Despite that, the task 
is not simple, since we need to minimize a functional we also need the gradients of 
this functional with respect to the variational parameters, i.e. the matrix elements 
appearing in the trial Hamiltonian Hq. The true minimum for F is obtained when 
Hq is the exact many-body Hamiltonian and in this case F[H] = —^Tr{Pe~^^). 

As an alternative approach we can approximate the propagator e~^^ appear- 
ing in the free-energy minimum, with some propagator that at the mean-field level 
reproduces the HFB approximation. We can then make use of the approach pro- 
posed almost 30 years ago in ref. [8], based on a formulation of the functional 
integral expression of the partition function different from the one obtained using 
the Hubbard-Stratonovich transformation. As well known, functional integrals 
for the partition function obtained by approximating infinitesimal evolution op- 
erators e"*^^ using the Hubbard-Stratonovich transformation do not lead, in the 
mean-field approximation, to the Hartree-Fock or to the Hartree-Fock-Bogoliubov 
approximation. In ref. [8] it has been suggested that adding appropriate terms in 
the elementary propagators appearing inside the functional integral, that do not 
change the functional integral, but that crucially change the mean-field, one can 
obtain either the Hartree-Fock or the Hartree-Fock-Bogoliubov approximation (cf . 
refs. [8], [9] for a detailed discussion). In ref. [8] only the unprojected partition 
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function has been discussed. However it enables to define which terms must be 
added inside the functional integral in order to obtain the HFB approximation by 
maximizing the integrand which contains only one trace (in contrast with several 
traces of the previous method). We shall rewrite here the functional that needs to 
be maximized. We use notations close to the ones in ref. [8]. i,j,k,l,.. denote 
single-particle states. The effective Hamiltonian is written as 

= < > alttj + 2 X/ ^ > alajdiak (6) 

ij ijkl 

K is the one-body term that also includes the chemical potential term —fip (for 
simplicity we do not distinguish neutrons and protons) and V is the two-body 
interaction. Let us define the operators 

Pij = a\aj, fjlj = a\a], rjij = ajUi (7) 

and the following matrix multiplication rules 

avp = X (^ijVijkiPki (8) 

ijkl 

l_iuf)'^ = J2 PijUijkifjii (9) 

ijkl 

fwl) = ^VijWijkif]ki (10) 

ijkl 

then the THFB approximation is obtained by maximizing e"^^ where 

12(4 + 4 + 4) -^^^ Tt{u) + ^iFA (11) 

where 

U = Q-^'^ (12) 

and 

% — < i\K\j > alttj — avp — {/j, + iv)wff — {/j, — iv)ui) (13) 

ij 
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Note here that the trace operation does not include the projectors. Also we have 
replaced u* in ref. [8] with w for more generality. In order to obtain the HFB 
approximation (without the projectors) one has to ^xv,u,w such that 

{vv)ijki = -V,lj^i, {wu)ijki = ~V^M (14) 

(note the peculiar order of the indices in the first equation) where is the an- 
tisimmetrized interaction. This is the only result we need from ref. [8]. A simple 
way to determine the matrix v is the following. Define first the indices a — {i, j) 
and 13 = {k, I) and the real hermitian matrix Wa,(3 = Vj^ji. Then decomposing 
W into its eigenmodes we have W = XuX where u is diagonal. Therefore one 
can make the choice v — \/—uX. Similarly we can decompose (with this 

grouping of the indices) in its eigenmodes V"^ = YaY with a being diagonal, and 
take u = w = ^—^aY. Since the eigenvalues in both cases can have any sign 
u, V, w are generally complex and therefore the Hamiltonian H of eq.(13) is not 
necessarily hermitian, except at the stationary point. 

It is sensible to consider the following projected functional 

^P=\ + 4 + 4) - I InTriPu) + i^fA (15) 

where U is given by eqs.(12)-(13). Once this functional has been minimized with 
respect to aij,jj,ij, Vij one can evaluate whatever thermal averages one wishes to 
compute. Hence the central task is to evaluate as efficiently as possible In Tr{PU) 
and its gradients with respect to Oij, jiij, Vij. Once this is done we can minimize 
f2p using for instance quasi-newtonian minimization techniques (ref. [10]). In 
the following section we shall discuss precisely how can one evaluate efficiently 
In Tr{PU) and its gradients. It will be shown that actually one can work directly 
in a canonical ensemble rather than in the grand-canonical one. As a final com- 
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ment, note that the free energy obtained in this way may not be optimal in the 
sense of the first method (this is the reason why we have used a different symbol), 
but at least we avoid entirely the problem of evaluating the logarithm of a pro- 
jected operator. Also, since is complex, it is best to minimize its real part and 
imposing a penalty function to suppress its imaginary part. That is it is better to 
minimize the functional G = Re{Vtp) + c{Im{Vtp))'^, c being a positive number 
that controls the strength of the penalty function. 

The strategy of approximating the exact propagator with an approximate one, 
and then evaluating its associated free energy, is not a new one. Also ref. [5] 
where the propagator was replaced by the projected SPA propagator, follows the 
same line, instead of minimizing the functional of eq. (1). 



3 Technique for the evaluation of In Tr ( P its gra- 
dients and expectation values in a canonical en- 
semble. 

Let us first rewrite the propagator U — exp{—/3Ti) in a more convenient form as 

U = exp(C) exp(^7,S7c) = exp(C)W (16) 

where 7c = col{a, a^) is the column vector of all creation and all annihilation 
operators for all Ng single-particle states (including the isospin label), and % = 
row(a^, a). The matrix 5" is a 2 AT, x 2Ns matrix of the type 




(17) 



with 

Su^-p{K-av), Si2^2p{ii + iiy)w, S21 ^ -2p{ii - iu)u (18) 

and 5*22 = — •S'n. The matrices 5*21 and 5*12 are skew-symmetric. Moreover C — 
J2i{Sn)ii. The most generic projector has the following structure 

P = ^L>*(a;)e'^^~ (19) 

where a; is a set of rotation angles (for example the Euler angles and/ or the fugaci- 
ties) and D* ia function of the rotation angles (for example the Wigner functions in 
the case of angular momentum projection) and r is a matrix in the single-particle 
indices. Therefore we will be concerned to the evaluation of quantities of the 
type Tr(e'^^™ W). In the following subsections we first recall briefly some of 
the properties of products of operators like W, called statistical density operator 
(SDO) which are not necessarily hermitian, and then we will apply them to the 
evaluation of the traces and their gradient with respect to the matrix elements of 
S 

3.1 Basic properties of Statistical Density Operators. 

For a thorough discussion we refer the reader to ref.[ll]. The most useful 
properties of the SDO's are the following. 

1. Group property. The product of two SDO's is a SDO with some matrix S 
associated to it. 

2. To any SDO there is a matrix associated to it. More precisely if 

W = exp(^7rS7c) (20) 

then the associated matrix is 

W^e^ (21) 
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(without the caret), which preserves the group structure, that is, if W — W1W2 
then W = W1W2. 

(0 l\ 

3. The matrix W satisfies the identity W ^ — rjWrj where rj — 

4. The matrix W can be diagonalized by a similarity transformation of the same 
type T, that satisfies T'^ = rjfrj (cf. ref. [12] for a proof). 

5. The eigenvalues of the matrix S come in opposite pairs (Ai, A2, ..Ajv^, — Ai, — A2, 
Xns)- Moreover the Grand-Canonical trace of W is given by (cf.ref.[12]) 

TrgcW = n (e^^/^ + e-^'/^) (22) 
fe=i 

6. Taking the square of eq.(22) we have 

(TrgcWy = det(l + W) (23) 

7. Any SDO can be written as a product of three special SDO's (ref.[l 1]), i.e. 

W = exp(^atCat)exp(^7, |^ ^ j 7e)exp(iaDa) (24) 

with 

C^WuW22\ D^W22^W2i, e-^=W,, (25) 

W21, W22 are the Ng x Ng blocks of W. Because of the relation Wr]Wr] = 1, 
the matrices C and D are skew-symmetric. 

The following properties are more relevant to the topic under discussion and 

can be proved using the above. 

8. The vacuum expectation values of W is given by det(H^22)^^^ and the vacuum 
expectation value of the product l^(l)iy(2) is given by 

< 0\W{1)W{2)\0 >= det{W22{l)y/^det{W22{2)y/^l['{l + i^k) (26) 



where the vj. are the eigenvalues of £)(1)C(2) which come in degenerate pairs and 

the ' denotes one eigenvalue per degenerate pair. 

9. The grand canonical trace of e"^ W can be rewritten using the eigenvalues jik 
of the following matrix 

( W^2 C \ 

M = 

V -D ) 
which come in degenerate pairs (/xi, /X2, ••, A*iv„, A*i, A*2, ••, A*ivJ, as 
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Tigele^^'Pr) = det(Pr22)'/' (28) 

k=l 

where z = e°, and the product includes one eigenvalue per degenerate pair. This 
very important property has been proved in a slightly different form in ref. [12]. 
The only difference here is the use of M = SpS~^, instead of M = S~^Sp (cf. 
ref. [12] for the details and the definitions of these matrices). In the next subsection 
we will discuss how the above equation allows us to work directly in a canonical 
ensemble rather than in the grand-canonical one. 

3.2 Evaluation of traces. 

Let us discuss first some of the implications of eq. [28]. Some of these have 
been discussed in ref. [12]. First if we know the matrix S that defines the SDO we 
can unambiguously determine the phase of det(H^22)^^^- In fact we can set 2; = 1 
in eq. (28) (that is a = 0) and obtain unambiguously det(1^22)^''^ using eq. (22) 
and eq. (28). 

Second, once we have done this, we can isolate the coefficient of in the prod- 
uct appearing in eq.(28), ^jv, and therefore the canonical ensemble trace for N 
particles of W is given by 
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The expression for the coefficient is the following 

Civ= Yl l^iil^i2-l^iN (30) 

ii<i2<---<iN 

These coefficients can be constructed iteratively by defining the coefficients ^ (n, s) 
for n particles using the first s distinct eigenvalues 

s) = C{n, s-l)+ i^sCin - 1, s - 1) (31) 

There is a nontrivial consequence of eq. (29). Qualitatively, the M matrix is 
scale independent. In order to see this, let us recall that the eigenvalues A are 
proportional to the temperature. Let T be the matrix of the eigenvectors of S that 
is 

f \ n \ 

(32) 

(cf. property 5.). We can order the eigenvalues so that the first Ng have positive 
real part. Let us call and la, 13 the partitions of T and (a, (3 = 1, 2). Also 
set D = and d = which represent the large and the small scale. Then we 
have 

VF22 = T21DI12 + ^22(^/22 = -^22 + >S'22 (33) 

W12 = TnDIu + Ti2dl22 = Lu + S12 (34) 

W21 = TnDhi + T22dh2 = L21 + ^21 (35) 

which shows the large and small scale ( arrays L and S) of the various blocks of 
W. One can write 

D = {1 + L^^S22)-^L^^L2^{1 + L-^S2i) (36) 

Using the definitions in eqs.(33)-(35) we have L22L21 — Ivihi- therefore, espe- 
cially at low temperature, D is scale (i.e. temperature) independent. Similarly for 
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C we have 

C = (1 + 5'i2iji2^)-^12-^22^(-'- '^'22-^22^) ^ (37) 

Again we have L12L22 — T11T21 which is temperature independent. The re- 
maining factors in C and D have a mild temperature dependence. Therefore the 
matrix M is roughly scale independent. Hence we arrive at the conclusion that 
the temperature dependence of the canonical ensemble trace of W is mostly in the 
vacuum expectation values ( the ratio in eq. (29) ). We obtained this conclusion 
under the assumption that the inverses of T^^^ and 1^,1} exist. This is not true in 
general. That is, the various block can be singular even though the T and / matri- 
ces are not. Despite this limitation in the proof, it is a surprising result. This scale 
independence of the the eigenvalues of M has been studied in detail in ref. [13] 
although only for the pure pairing model. 

This nearly temperature independence implies that at low temperature does 
not contribute to the energy at low temperature. This in turn implies that any 
dependence of the energy from the number of particles must come from the vac- 
uum contribution. Differently stated, the chemical potential is the only parameter 
that determines the dependence of the energy on the number of particles. The 
chemical potential, which plays an essential role in the grand-canonical ensem- 
ble, retains its importance also in the canonical ensemble. Being a free parameter 
in the canonical ensemble, and since we seek to minimize the free energy func- 
tional, it must be fixed in order to minimize the free energy, hence d^jQ. = with 
Vt given by eq. (15). We have assumed that the eigenvalues of M are dominated 
by the C and D blocks. If, as in the case of lack of pairing, these matrices are 
zero, the matrix 1^22^ which appears in M will be relevant in the determination of 
the energy and this matrix has an exponential dependence on the temperature. 
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So far we have discussed the evaluation of the trace in the canonical ensemble, 

without any other projector. Let us turn now to the evaluation of Trc(exp(a^ra)W). 
Note that < 0|exp(a^ra) |0 >= 1. Let us first rewrite 

1 1|^"0\ 1 

exp(a^ra) = exp(- ^rii)exp(-7, I I 7c) = exp(- ^rii)R (38) 

and make use of property 8. The result for the vacuum contribution is det(iy22)^^^ 
as in the case r — 0. The associated matrix to 



W'^RW 



(39) 



IS 



, e'" 

w ^\ \w 

e-^~ 



(40) 



evaluating the matrix M' of W (cf. eq.[27]) we obtain that the eigenvalues of M' 
are the same of the matrix 



e'' ^ 



V 



M 



(41) 



e*^ 

where M is relative to W only. The final result for the canonical trace in this case 
is 

Tre(exp(atra)W) = det(W22)'/'eN(/^^'^) (42) 

where the //^^^ are the eigenvalues of M^^^ and ^iv(A*^''^) is given by the recursion 
relation of eq.(31) with the "rotated" eigenvalues iJ,^'''K The vacuum contribution 
det(H^22)^'^^ is again given by 



det(W^; 



22 



1 1/2 



n£i(e^'/^+e-^'/^) 



(43) 
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Where the //^'s are the unrotated eigenvalues of M. We stress again that the 
eigenvalues of M'^^^ are two-fold degenerate and that only one eigenvalue per pair 
must be taken in the evaluation of ^n- 

From a numerical point of view, the determination of the eigenvalues ( A^, — A^) 
of S (of. eq.(32)) does not pose a problem. The determination of the eigenvalues 
il^k, l^k) and {iJ,'^^\ iJ,^^^) does. As discussed previously, the partitions of W carry 
both a large and a small scale that cancel out in the determination of the arrays 
C and D. This analysis however is based on the existence of the inverses of the 
partitions of T and of its inverse, which cannot be guaranteed. The most reliable 
way to prevent loss of accuracy in the determination of M for arbitrarily large 
values of /3 is the following. Let us divide the interval [0, /3] in Np equal intervals 
and let us define the matrix u associated to the propagator in each interval. The c 
and d (we use small letters for each interval) matrices can be determined without 
loss of accuracy. In force of the group property of the propagators we have (right 
to left propagation) 

W{n) = uW{n - 1) (44) 

where n — 2,..Np and W{1) — u. We seek a "propagation" law for C and 
D as the interval is enlarged. Define the following auxiliary matrices for each 
elementary interval 

7 = Tin 1*12, 5 = ^i2l^in (45) 

Then we have 

C(n) = Mil [C(n - 1) + 7][1 + dC{n - l)]-^u:^^ (46) 

At each step, C{n) is well determined since the matrix u is close to 1. For the 
determination of D we found more convenient to propagate W from the left to 
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right, i.e. 

W'{n) = W'{n + l)u, {n = N -1,N -2,..l) (47) 
with W'{Nj3) — u, using the definition of D we find 

D{n) = -^22^1 + D{n + l)c]-^[D{n + 1) + S]uu (48) 

with D{Np) — d. Again in the backward propagation, D is numerically stable. 
There is a final point which must be discussed, that is whether U22 has an in- 
verse. If we were working with hermitian S the Bloch-Messiah theorem states 
that indeed, 1^22 can have zero eigenvalues. These however would be associated 
to empty single-particle states, which are unlikely to happen, since the preferred 
basis is the harmonic oscillator basis. 

The remaining task is to give some prescriptions for the determination of the 
gradients of Trc(RW) with respect to all matrix elements Sij. It is simpler to 
evaluate the variations of ln(Trc(RW)), and the result is 

5(ln(Tre(RW))) = ^ T^^,^^^^^) ' E TTTT^^'"'^) + E fk(/x«)5(/x?) 
j^^^ i + e *= 2 k=i 1 + A^k k=i 

(49) 

where the fractional partition function (fpf) 

= ^?(!« (50) 

satisfy the sum rule Z^^i Hkfk — N (the number of particles). These fpf fk can 
be obtained from the by excluding the eigenvalue Hk and decreasing the number 
of particles by one. The various ^A^, Sfik etc. can be evaluated using first order 
perturbation theory. Perturbation theory on M or M^'^^ is a bit involved. The 
simplest way is first to determine 5Wij at each interval . To simplify the notations 
let us call with greek letters the eigenvalue indices of S and let T be the matrix 
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that diagonalizes S. Using the sum convention over repeated indices, we have 

5Wij = TiJW^^T^l (51) 

SWa^ = Fa/jSSa/S nO SUHl (52) 

SSap = T^^^SSrsTsp (53) 

where 

= — , (Xa ^ A^) (54) 

and = e'^". The matrix bS^s has only two non-zero elements. We have 
performed a numerical comparison of this method to evaluate the gradients with a 
simple-minded approach that uses the numerical evaluation of the derivatives and 
found that it is numerically stable even for large values of /5 for not too large iVg. 

3.3 Evaluation of one-body and two-body expectation values. 

Let us now turn to the problem of evaluating canonical ensemble traces in- 
volving one-body and two-body observables. Consider the operator 

Q{e) — exp(ea^qa) (55) 

where a)qa — J2i o^Qij^^j' let us evaluate 

Y,{e) = Tre(QRW) = Yq + eYi + ^e^Ys + ... (56) 

where the subscript c stands for canonical, up to second order in e. The second 
order term gives the matrix expectation values of (a^ga)^. a^qa could be for in- 
stance one of the eigenmodes of alaiv{ij, kVja^^ak in terms of the pairs (i, I) and 
(j, k). Let R be the "rotation operator" that appears in the remaining projectors. 
Consider 

Yg,{z, e) = TVge(e"^ QR W) = Yg,o + e Y,,i + ^e' Y,,^ (57) 
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where z — e". We simply have to deterniine the coefficient of in YgJyZ^ e) up 
to second order in e for particles. From eq.(42), we have 



(58) 



i=l 



where the product contains the eigenvalues (one per degenerate pair) of 

„ i QR ^ \ 
M = M 

V RQ ] 

where R — e^. The matrix M" can be rewritten as 

/ 1 n \ ( r^ n \ / 



(59) 



M = 



1 
R 



1 



R* 



(60) 



Q 

V g 

using the fact that it! is unitary in the cases of physical interest. In order to simplify 
the notations let us set 



(61) 



and 



r/ = (1 + zM('-))-^5M('-) (61) 
Then, expanding the determinant 1 1 + zM" \ up to second order in r] we obtain 

1 1 + zM" I = 1 1 + zM^'^ I (1 + ]-ii{ri) + itr^ (r;) - itr(r;2^ (62) 

2 8 4 

Using the diagonal representation of M*^**) and the fact that Ygc must be a polyno- 
mial in z we obtain, after some algebra 

y^c. = r22r/^ E f 11(1 + ))) ^A^r)[(5o.. + (63) 

3 \«7^J / 



and 



i;c2=r22r/^[E n (i+^4^^)hvrvfG,,+ 
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E f 11(1 + zf,f^)\ zf.P + {d,)jj] (64) 

where 

(65) 

In eqs. (63),(64), sum and product are performed considering one member per 
degenerate pair and i represents the state which has the same eigenvalue of the 
state labeled by i. Since terms containing 1/(1 + z[i(^)'^ must cancel out, we find 
(5i)m = (^i)ii aiid (5i)j| = 0. We can now proceed directly to the canonical en- 
semble and define in addition to the fractional partition function /j of the previous 
subsection, the fractional partition function for a pair {%]) 

fS^ - e^;^^^-^^^^'^ ^ (/^r, /^f )) (66) 

Appropriately i^f"^ fi and i^f^ (Ji'p flf can be interpreted as the occupation numbers 
and pair occupation numbers of nucleons.The quantities i.i;[\i'p fl"^^ satisfy the 
sum rule Y,i<j fJ't'' l^'t'' fif ~ N{N — l)/2. The final results for one-body and 
two-body canonical ensemble traces are as follows 

IVe(atqaRW) = tr,(RW) ^/.«fk(//«)((52)jj + (62^) (67) 

k 

and 

IV,((atqa)2RW) = tr,{RW)[^ t^^^HS^h + {32^] + E /^i/^j^if Gy] (68) 

i i<j 

The first term can be eliminated by rewriting (a^ga)^ = a^q^a + alajquqj^kO^kai, 
hence 

Tre(a|a|qiiqjkakaiRW) = tr,(RW) ^ /x^Vf ^Gy (69) 

i<j 
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By extending the definition of as fij = /^j = fjj = fjj we obtain 



Tre(ata|qi,qjkakaiRW) =tre(RW)^//S'-Vf ff [;^((^i)ii(<^i)jj - (<^i)ij(<^i)ji] 

(70) 

where now the sum is unrestricted, that is, it extends to all eigenvectors. The 
first term has the form of a direct contribution and the second of an exchange 
contribution. As a final point, the matrix elements are of the type 



-1 



where 



/ 



U 



% ^ 
[O q) 

\ 



u] 



^3 



1 

where V is the matrix that diagonalizes M'^'^\ i,e, 

//(^) 



V 



M^'^ = V 
the //('■^ being the eigenvalues. 



/x^*") 



(71) 



(72) 



(73) 



4 Conclusions. 

In this work we have discussed a variant of the thermal Hartree-Fock-Bogoliubov 
approximation (with the restoration of quantum numbers) whereby an approxi- 
mation to the exact free-energy is considered for a variational calculation which 
in the case of absence of the projectors gives the standard THFB. This variant is 
given directly in the canonical ensemble and equations for evaluation of traces, 
gradients of traces and one and two-body expectation values are given. 
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